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ABSTRACT 

Aims. We study the exponential scale length of total radio emission, the spectral index distribution, and the linear 
radio polarization in the Scd galaxy M33. 

Methods. Observations were carried out using the 3.6 cm dual channel and the 6.2 cm four channel receivers of 
the 100-m Effelsberg telescope along with the L-band VLA D-array at 20 cm. 

Results. High spatial resolution and sensitivity in both total and linearly polarized radio continuum emission from 
M33 were achieved. We found considerable extended emission, not only from the main arms IS and IN, but also 
from the weaker arms. The large-scale magnetic field exhibits well-ordered spiral structure with almost the same 
orientation as that of the optical spiral arms, however, it does not show a clear structural correlation or anti- 
correlation with the optical arms. There is a north-south asymmetry in polarization that is frequency-dependent. 
We found that the ring mean spectral index versus radius increases faster beyond _R = 4kpc. At each wavelength, 
the exponential scale length is larger inside than outside ii = 4kpc. 

Conclusions. From the larger scales lengths at i? < 4 kpc, we conclude that star forming regions are mainly spread 
over the region i?<4kpc without a dominant nuclear concentration. Furthermore, at i?<4kpc, a spatial corre- 
lation between cosmic rays and star forming regions may exist. From the behaviour of the mean spectral indices 
obtained from different pairs of the radio continuum data at 3.6, 6.2, and 20 cm, we confirm that a decrease in the 
thermal fraction causes an increase in the spectral index. The frequency-dependent asymmetry in the polarization 
hints to an asymmetry in Faraday depolarization. 
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1. Introduction 

M33(NGC598), the nearest Scd galaxy at a dis- 



tance of 840 kpc (1" ~ 4pc, iFreedman et al 
199ll ). has been extensively studied at radio 
wavelengths. Terzian fc Pankonin ( 19721 ) detected 
some diffuse emission at 318 and 606 MHz us- 
ing a ngular resolutions of 1 7' an d 10', respec- 
tively. Ilsrael fc van der Kruid (|l974l ) observed M33 
at 1415 MHz with the Synthesis telescope in 
Westerbork (WSRT). However, they did not find 
any extended spiral arm emission. This was first 
revealed by observations with the 100-m Effelsberg 
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dish at 6.2 cm (jvon Kap-Herr et al. _[l97i)_ that 
was extensively studied by BerkhuiisenI ( 1983 ). 
Radio pol arizati o n obs ervations of M33 were first 
made by lBeck.l (Il979l) a t 21.1cm and 11.1cm. 



Buczilowski fc Beckl ( 1987 ) observed this galaxy at 
17.4cm, 6.3 cm, and 2.8 cm. They detected polar- 
ization mostly in the northern half of M33 and 
proposed a bisymmetric magnetic fiel d structure 
The n ext interferometer observations ([Puric et al 



1993t IViallefond etldl Il986l: iGordon et all Il999^ 



were carried out in total intensity to study point 
sources like Hll regions and supernova remnants 
(SNRs) within M33. 

Generally, single-dish observations should be 
carried out to study the extended radio emission. 
So far, single-dish observations have been limited 
by sensitivity because M33 is relatively weak in ra- 
dio emission, especially in polarized emission (with 
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Table 1. Positional data on M33. 

Nucleus position (J2000)' RA = l''33'"51.0" 

DEC = 30°39'37.0" 
Position angle of major axis^ 23° 
Inclination^ 56° 
Distance"^ (l"=4kpc) 840 kpc 

^ de Vaucouleurs fc Leach (1981) 
^ Deul fc van der Hulst f 1987) 
^ Regan fc Vogel ('19941 
Freedman et al. (1991) 



Table 2. Characteristics and performances of 
3.6 cm and 6.2 cm receivers of the 100-m Effelsberg 
telescope. 



an integrated flux density of '^83mJy at 6.2 cm, 
Buczilowski fc Beck Ill99l1) and because of M33's 
large extent in the sky (~ 35' x 40'). 

This paper is the first of a series providing a de- 
tailed investigation of the radio continuum emission 
of M33. We describe new observations of the radio 
continuum emission at 3.6 cm, 6.2 cm, and 20 cm 
with much improved detection limits of total and 
polarized intensity and with higher spatial resolu- 
tions. These data allow a much more detailed study 
of polarization and magnetic field within M33 than 
was possible before. We discuss the distribution and 
degree of linear polarization, integrated fiux densi- 
ties and distribution of the total spectral index (ob- 
tained from total intensities of the radio continuum 
emission including both thermal and non-thermal 
emission). We also compare the radial, exponential 
scale lengths of the radio emission at different wave- 
lengths. 

A comprehensive treatment of the separation of 
thermal and non-thermal emission, rotation mea- 
sures, Faraday effects, depolarization, magnetic 
field strength and structure, radio-IR correlation, 
and cosmic ray distribution will be given in further 
papers of this series. 

The observations and data reduction are de- 
scribed in Sect. 2. Results including distributions 
of total intensity and polarization, integrated flux 
densities and total spectral index distribution are 
discussed in Sect. 3. Discussion and concluding re- 
marks are presented in Sect. 4. 

2. Observations and data reduction 

We performed both single-dish and interferome- 
ter observations. The single-dish observations at 
3.6 cm and 6.2 cm were made with the 100-m 
Effelsberg telescop^. The interferometer observa- 

^ The 100-m telescope at Effelsberg is operated by 
the Max-Planck-Institut fur Radioastronomie (MPIfR) 
on behalf on the Max-Planck-Gesellschaft. 
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83.6 


146 


Tsys (K) 


23-25 


25-28 


n/S. (K/Jy) 


2.62 


2.4 



tions at 20 cm were carried out with the Very Large 
Array (VLA|). 



2.1. Effelsberg single-dish observations 

We observed M33 at 3.6 cm during several periods 
from August 2005 to March 2006 with the broad- 
band system (bandwidth ~ 1100 MHz) centered at 
8350 MHz. The receiver is a dual-channel correla- 
tion radiometer with cooled HEMT ampliflers and 
a cooled polarization transducer located in the sec- 
ondary focus of the 100-m telescope. The system 
is set up for receiving Left Hand Circular (LHC) 
and Right Hand Circular (RLC) polarization. M33 
was scanned along RA and DEC with a scanning 
velocity of 50'/min on a grid of 30". The scanned 
area of 40' x 50' was centered on the nucleus of 
the galaxy. The total observing time was about 100 
hours. We obtained 56 coverages with a typical sen- 
sitivity per scan of 2mJy/beam area in total inten- 
sity and 0.4mJy/beam area in polarization. The 
Hah Power Beam Width (HPBW) at 3.6 cm is 83.6" 
which corresponds to a linear resolution of about 
335 pc. 

The sources 3C 48, 3C 286, and 3C 138 were used 
for pointing and focusing during the observations. 
Calibration of the flux density and polarization an- 
gle was achieved by observ ing 3C286. We use d the 
flux density scale given bv lBaars et al.l (jl977l ). 

The data reduction was performed in the N0D2 
data reduction system (|Haslamlll974l ). In order to 
remove the scanning effects due to ground radia- 
tion, weather condition and receiver instabilities, 
we applied the sca r ining removal program, Presse, 
of Sofue fc Reich (|l979l) . The r.m.s. noise after 
comb ination of the coverages (jEmerson fc Graeve 
Il988l ) is -220AtJy/beam in Stokes I (total inten- 
sity) map and ~70/iJy/beam in Stokes U and Q 
(polarization) maps. The final maps in U and Q 



The VLA is a facility of the National Radio 
Astronomy Observatory. The NRAO is operated by 
Associated Universities, Inc., under contract with the 
National Science Foundation. 
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Fig. 1. Total intensity 
emission from M33 at 
3.6 cm observed with 
the Effelsberg telescope, 
with apparent B-vectors 
(E-vectors rotated by 
90°) of polarized in- 
tensity superimposed. 
The HPBW is 120" (the 
beam area is shown in 
the left-hand corner). 
Contour levels are 0.5, 
1, 2, 3, 4, 8, 16, 32, 
64mJy/beam. The rms 
noise is 0.14mJy/beam 
in total intensity and 
0.04mJy/beam in linear 
polarization. The vectors 
were plotted starting 
from 0.1 mjy /beam 
(2.5 a). A vector length 
of 1' represents a po- 
larized intensity of 
0.5 mjy/beam. 



were combined to produce maps in polarized in- 
tensity and in polarization angle, correcting for the 
positive bias in po larized intensity due to noise 
( Killeen et al. 198i). ?igs. 1 and 2 show the re- 
sultant maps of the total and polarized intensities, 
respectively, smoothed to an angular resolution of 
120". 

Observations at 6.2 cm were carried out in the 
summer of 1996 using the 4850 MHz HEMT re- 
ceiver installed in the secondary focus of the 100-m 
Effelsberg telescope. This is a 4-channel 2-beam sys- 
tem with cooled HEMT pre-amplifiers (see Table 
2 for more information about this receiver). M33 
was scanned along its major axis and perpendicu- 
lar to it on a grid of 1' with a scanning velocity 
of 907min. The field size is 100' x 80'. The num- 
ber of useful coverages is 12 for Stokes I maps and 



14 for each U and Q maps. The 6.2 cm data reduc- 
tion was performed similar to that at 3.6 cm. The 
r.m.s. noise after combination is ^ 700 /iJy/beam in 
Stokes I (total intensity) map and 180/iJy/beam 
in Stokes U and Q (polarization) maps. The final 
maps of total and polarized intensities (smoothed 
to 180") are shown in Figs. 5 and 6, respectively. 

2.2. VLA observations 

We observed M33 at 20 cm (L-band) in two 50 MHz 
bands centered at 1385 and 1465 MHz with dual 
circular polarization in the continuum mode of the 
VLA D-array. The observation dates were 06-11-05, 
07-11-05, 08-11-05, 13-11-05, and 06-01-06. To cover 
the central 30' x 45' with equal sensitivity, we made 
a mosaic with 12 pointings with a spacing of half 
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Fig. 2. Linearly polarized emission (contours and grey scale) and degree of polarization (vector lengths) 
from M33 at 3.6 cm observed with the Effelsberg telescope, smoothed to an angular resolution of 120" 
(the beam area is shown in the left-hand corner). Position angles of the vectors show the orientation of 
the apparent magnetic field component perpendicular to the hne of sight. Contour levels are 0.1, 0.2, 0.3, 
0.4mJy/beam. The rms noise is 0.04mJy/beam. A vector length of 1' represents a degree of polarization 
of 12.6%. 



the primary beam width, 15', (3 pointings along RA 
and 4 pointings along DEC). The coordinates of the 
12 pointings are listed in Table 2. The observations 
were made in cycles with 6 min on each pointing of 
the mosaic. This observing cycle was repeated 17 
times. At the beginning and end of each 4-cycle, 
we observed 3C 138 and 3C 48 as the primary flux 
density calibrators. The antenna gains and phases 
were calibrated every 43 min with the phase cali- 
brator 0029 + 349. The source 3C 138 was also used 
for polarization calibration. 

The standard procedures of the Astronomical 
Image Processing System (AIPS) were used to re- 
duce the data. The data from each day were cali- 
brated separately to ensure that there were no day- 



to day amplitude discrepancies. Then they were 
split into 12 single-source visibility data sets corre- 
sponding to the different pointings. After flagging, 
we combined the source visibility data from differ- 
ent days of the same source. We self-calibrated the 
combined visibility data in both phase and ampli- 
tude & phase. Final images were made with IMAGR 
parameter Robust=0. The rms noise of different 
pointings in both total intensity (I) and polarization 
(U and Q) arc listed in Tabic 3. We convolved all 
the single images to the largest clean beam of 51.4" 
in both RA and DEC, the individual images were 
reprojcctcd on a single grid centered at RA = Ol'' 
33" 51.41^ and DEC = -F30° 46' 45.00". Then, we 
combined the single images to a mosaic with equal 
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Fig. 3. A sketch of the optical arms ( Sandage fc Humphrevsl 198(]f ) overlayed on the total (left panel) and 
linearly polarized (right panel) intensities at 3.6 cm. Grey scales and Contour levels are the same as in 
Figs. 1 and 2. 



Table 3. Central positions of the 12 fields observed with VLA at 20 cm. 



Pointing 


RA(q) 


DEC (5) 




I rms noise 


U (Q) rms noise 


HPBW 


# 






C ' ") 






(/xjy/beam) 


(/^Jy/beam) 


(arcsec^) 


1 


01 


35 01.26 


+30 16 


45 


00 


26 


23 


50.7 X 49.5 


2 


01 


33 51.79 


+30 16 


45 


00 


23 


20 


50.2 X 46.1 


3 


01 


32 42.31 


+30 16 


45 


00 


27 


18 


50.6 X 43.4 


4 


01 


35 01.26 


+30 31 


45 


00 


26 


21 


48.1 X 44.0 


5 


01 


33 51.60 


+30 31 


45 


00 


28 


20 


47.8 X 45.4 


6 


01 


32 41.96 


+30 31 


45 


00 


32 


19 


51.4 X 42.2 


7 


01 


35 01.26 


+30 46 


45 


00 


26 


21 


50.0 X 46.1 


8 


01 


33 51.41 


+30 46 


45 


00 


34 


21 


49.2 X 42.4 


9 


01 


32 41.58 


+30 46 


45 


00 


23 


19 


51.1 X 44.0 


10 


01 


35 01.26 


+31 01 


45 


00 


36 


23 


48.5 X 45.2 


11 


01 


33 51.21 


+31 01 


45 


00 


24 


21 


51.1 X 43.3 


12 


01 


32 41.21 


+31 01 


45 


00 


21 


19 


48.0 X 44.0 



weights using the task 'LTESS' which also corrects 
for primary beam attenuation. 

There is a considerable lack of flux density in the 
VLA maps: the integrated total flux is only ^ 15% 
in total intensity and ^ 55% in polarized intensity 
of those in the Effelsberg maps (determined as de- 



scribed in Sect. 3.3). This is not surprising because 
of the missing short spacings in the forms of severe 
attenuation of angular scales larger than ~ 13' (cor- 
responding to the shortest spacing of ~ 40 m) and 
the negative bowl artifact around strong sources. 
To correct for the missing spacing, we combined 
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Fig. 4. Apparent B-vectors (E-vectors rotated by 90°) at 3.6 cm superimposed on an optical image (B- 
band, taken from the STScI Digitized Sky Survey) of M33. A vector length of 1' represents a polarized 
intensity of 0.3mJy/beam. The vectors were plotted starting from 1.5 tr. 



the VLA map with t he Effelsberg single-dish 20 cm 
map (|FletcherJl20"0lh in the uv plane. To avoid dis- 
tortions in the combined total intensity (I) map 
due to strong sources, we first subtracted these 
sources from the VLA map. The sources were iden- 
tified in the Effelsberg map by comparing with the 
VLA map smoothed to the Effelsberg resolution. 
Then they were also subtracted from the Effelsberg 
map. The residual VLA and Effelsberg maps were 
combined using the task TMERG' with a maxi- 
mum spacing of 0.236 kA for the Effelsberg map 
and a minimum spacing of 0.170kA for the VLA 
map in total intensity. The corresponding maxi- 
mum and minimum spacing adopted in polarizar- 



ion(U and Q) are 0.171 kA and O.lSOkA, respec- 
tively. These maximum and minimum spacing val- 
ues were chosen so that the integrated flux densities 
of the merged maps became equal to those of the 
Effelsberg maps. After combination, the VLA point 
sources were added again. The r.m.s. noise in the 
combined maps are ~70/iJy/beam area in Stokes 
I and ~ 25 /iJy/beam area in Stokes U and Q. The 
resultant total and polarized intensity images are 
shown in Figs. 7 and 8, respectively. 
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Fig. 5. Total intensity emission from M33 at 6.2 cm observed with the Effelsberg telescope, with an 
angular resolution of 180" (the beam area is shown in the left-hand corner). Apparent B- vectors of the 
polarized intensity (E-vectors rotated by 90°) are superimposed. Contour levels are 2.2, 4.4, 6.6, 8.8, 
17.6, 26.4, 35.2, 70.4mJy/beam. The rms noise is 0.56mJy/beam in total intensity and 0.12mJy/beam 
in polarization. The vectors were plotted starting from 2.5 a. A vector length of 1.5' represents a polarized 
intensity of ImJy/beam. 



3. Results 

3.1. Total intensity 

Our total intensity map of M33 at 3.6 cm is a 
fa ctor of ~ 6 more sen s itive than the 2.8 cm map 
of iBuczilowski fc Beckl (jl987() . The map at 3.6cm 



with the ori ginal angular reso l ution of 84" was 



presented by iTabatabaei et al.l (|2007l ). Even after 



smoothing to 120", the total intensity map (Figs. 
1 and 3) shows that not only the two main opti- 
cal arms I ^ and I N are well pronounced, but also 

^ We refer to the arm notation of 
iHumphrevs fc Sandagd (|l980l ): IN to VN for the 
northern and IS to VS for the southern arms. 



the northern arms UN, IIIN, IV N, and VN are 
(at least partly) visible at 3.6 cm. The total inten- 
sity of the extended emission is about 4mJy/beam 
in the arm IN which increases to 8mJy/beam in 
parts of the arm IS. The intensity is smaller than 
3mJy/beain in the other arms. The flux densities 
of 11 brightest HII regions at 3.6 cm are presented 
in Table 4. 



At 6.2 cm, the rms noise (0.7mJy/beam) is im- 
proved by a factor of two with respect to the pre - 
vious observations by IBuczilowski fc Beck ( 1987 ) 
The northern arms and I S plus parts of the south- 
ern arms II S and III S are visible in the second con 
tour in Fig. 5. 
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Fig. 6. Linearly polarized emission (contours and grey scale) and degree of polarization (vector lengths) 
from M33 at 6.2 cm observed with the Effelsberg 100-m dish, with an angular resolution of 180" (the beam 
area is shown in the left-hand corner). Position angles of the vectors show the oricntaiton of the magnetic 
field component perpendicular to the line of sight. Contour levels are 0.3, 0.6, 0.9, 1.2, 2.4mJy/beam area. 
The rms noise is 0.12 mJy/beam. A vector length of 1.5' represents a degree of polarization of 11.7%. 



With a linear resolution of ~ 200 pc 
(HPBW = 51.4") and a sensitivity of 70^Jy/beam, 
a large number of point sources are revealed in the 
total intensity map at 20 cm (Fig. 7). A list of the 
flux densities of the radio point sources at 20 cm 
at 7" resolution was presented by Gordon et al 



(|l999lFI . The contour levels shown in Fig. 7 starts 
at the 4(7 level. The strongest extended emission 
is found in the arm I S with a flux density of 
> 3.6 mJy /beam, followed by parts of the arm IN 
and the extended central region of the galaxy. At 
20 cm, the extended emission is detected up to 

They observed M33 with the VLA and WSRT at 
1.4 GHz, measured the flux density of point sources, but 
missed the extended emission. 



larger radii and is more radially symmetric than at 
3.6 cm and 6.2 cm. 



3.2. Distribution of polarized emission 

The polarized emission at 3.6 cm is shown in Fig. 2 
in grey scale and wi th contours. A sketch of t he op - 
tical spiral arms by Sandage fc Humphrevd ( 1980l ) 
superimposed on the polarized emission is shown 
in Fig. 3. Strong polarized intensities (more than 
0.4mJy/beam) are found in between and parts of 
the northern optical arms IV N and IN (see Fig. 
3) and also in central regions of the galaxy. It 
seems that the polarized emisson, as a whole, obeys 
a spiral pattern structure. Whether this structure 
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Fig. 7. The combined VLA + Effelsberg map of the total intensity emission from M33 at 20 cm, with 
apparent B-vcctors of polarized intensity (E-vcctors rotated by 90°) superimposed. The HPBW is 51" 
(the beam area is shown in the left-hand corner). Contour levels are 0.3, 0.6, 0.9, 1.2, 1.8, 2.4, 3.6, 4.8, 9.6, 
19.2mJy/beam. The rms noise is 0.07mJy/beam in total intensity and 0.025 mJy/beam in polarization. 
The vectors were plotted starting from 1.5 a. A vector length of 51" represents a polarized intensity of 
0.1 mJy/beam. 
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Fig. 8. The combined VLA + EfFelsberg map of the Unearly polarized emission (contours and grey scale) 
and degree of polarization (vector lengths) from M33 at 20 cm. The HPBW is 51" (the beam area is 
shown in the left-hand corner). Contour levels are 0.07, 0.14, 0.21, 0.28, 0.42, 0.56, 0.84mJy/beam. The 
rms noise is 0.015 mJy/beam. A vector length of 51" represents a degree of polarization of 8.7%. 



Table 4. The 11 brightest HII regions at 3.6 cm. 



Object 


RA(a) 


DEC (5) 


Flux density 




m 


(° ' ") 


(mJy/beam) 


NGC604 


1 34 32.9 


30 47 19.6 


53.7 ±0.4 


NGC595 


1 33 32.4 


30 41 50.0 


20.7 ±0.7 


IC133 


1 33 15.3 


30 53 19.7 


9.4±0.5 


B690 


1 34 06.5 


30 47 40.0 


5.3 ±0.3 


B61/62 


1 33 44.8 


30 44 50.0 


4.9 ±0.5 


IC132 


1 33 17.1 


30 56 51.8 


4.7 ±0.2 


IC131 


1 33 12.5 


30 45 21.5 


4.7 ±0.3 


NGC588 


1 32 44.7 


30 39 20.6 


4.5 ±0.3 


IC142 


1 33 56.7 


30 45 51.8 


4.1 ±0.2 


B691 


1 34 15.4 


30 52 21.8 


4.1 ±0.5 


NGC592 


1 33 12.6 


30 38 51.3 


4.1 ±0.3 



mimics the optical arm structure or not can be 
seen in Fig. 3, where a ske tch of optical arms 



( Sandage fc Humph reys! 1980l ) is superimposed on 
the polarized intensity map. The polarized emission 
is found partly in the optical arms and partly in be- 
tween the arms. In the other words, there is no clear 
correlation between the optical arms and polariza- 
tion arm-like structure. The vectors in Fig. 1 (and 
also 2) show the intrinsic orientation of the ordered 
magnetic field in the plane of the sky (PI^B^g^j^), 
as Faraday rotation is not significant at 3.6 cm. A 
comparison between the orientation of the ordered 
magnetic field and optical structure is shown in 
Fig. 4. Again, it is seen that the ordered magnetic 
field traces out some sort of spiral structure which 
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does not coincide well with the optical spiral arms, 
although there is a general similarity in the orienta- 
tions. The strongest arm like filament of the mag- 
netic field exists in between the northern optical 
arms IV N and VN, but also covers parts of these 
two arms. Generally, in some regions the ordered 
magnetic field alignes parts of the optical arms (like 
parts of the arms IN and VN), or is driven to the 
optical interarms, or locally distorted in orientation. 
The latter is mostly seen in the southern part of the 
galaxy. 

The vector lengths in Fig. 2 represent the degree 
of polarization, P. Apart from in the northern arm 
IV N, P is not only high where the polarized inten- 
sity is strong: the southern optical arms IS, II S also 
show high P (larger than 20%). The highest degree 
of polarization, ~ 40%, is found near NGC604 at 
RA = 01'' 34'" 49.67'' and DEC = -H30° 47' 19.49" 
(Fig. 3). However, it is less than 5% within the giant 
HII regions (eg. NGC604 and NGC595). This plus 
a weaker ordered magnetic field (and low P) in the 
central southern arm (IS, the strongest arm in ra- 
dio, IR and optical bands) indicates that high star 
forming activity disturbs locally the regular mag- 
netic field and causes the patchy distribution of the 
degree of polarization within the arms. 

At the angular resolution of 180" at 6.2 cm, the 
magnetic arm like filaments are smoother and not 
well distinguishable (Fig. 6). The strong polarized 
emission (more than 1.2mJy/beam) again occurs 
in parts of the arms IV N and IN, but in contrast 
to that at 3.6 cm, not in the central regions of the 
galaxy (Fig. 5). The high P at 6.2cm (> 20%) has 
almost the same distribution as at 3.6 cm. In the 
arm IV N, the maximum P is about 30%, similar to 
that at 3.6 cm. 



3.3. Integrated flux densities 

After subtraction of the background sources 
(14 sources whose 20 cm flux densities were 
> 5 mJy/beam), we integrated the flux densities in 
rings in the plane of the galaxy around the center 
out to a radius of 8.5 kpc (beyond the area shown 
in Figs. 1 to 8)(see Table 1). The integrated flux 
density in total intensity, St, and in linear polar- 
ization, Sp7, at different wavelengths are given in 
Table 5. The corresponding absolute total flux den- 
sity errors were obtained taking into account the 
systematic error due to the baselevel uncertainty. 

The average degree of polarization (P), which is 
the ratio of the integrated flux density in polarized 
intensity and the integrated flux density in total in- 
tensity, is obtained at different wavelengths (Table 
5). For a consistent comparison, we also obtained 
the average degree of polarization after smooth- 
ing the 3.6 cm and 20 cm maps to the beam width 
of 180", P(180"). At 3.6cm and 20cm, P(180") is 
smaller than P due to beam depolarization. P(180") 
is almost the same at 3.6 cm and 6.2 cm, indicating 
negligible global Faraday depolarization between 
these two wavelengths. 

The integrated total flux densiti e s in to- 
tal intensity obtained by iBuczilowskil (|l988l ) of 
1100±167mJy at 6.3cm and 2990±440mJy at 
20 cm out to a radius of 12 kpc agree with our val- 
ues of 1300 ± 140 mJy at 6.3 cm and 3164 ± 88 mJy 
at 20 cm within the errors. Our obtained value 
at 6.2 cm i s even closer to that obtained by 
BerkhuiisenI (| 19831) (1300 ± 200 mJy). The inte- 
grated total flux density in total intensity found by 
Buczilowski (1988) at 2.8 cm is significantly smaller 
than that we found at 3.6 cm because the galaxy 
size is much larger than the ma ximum spacing 
of th e 2.8cm multibeam system ( Emerson et al.l 
1979h . 



Fig. 7 shows, at a linear resolution of ~200pc, 
how the apparent magnetic field vectors (E-|-90°) 
at 20 cm are distributed in the plane of the sky. At 
this wavelength, the vectors are strongly Faraday 
rotated and therefore are not oriented along the 
optical arms. The polarized signal is strongest in 
the northern part of the galaxy, where the order- 
ing in a spiral like pattern is very clear (Fig. 8), 
while the lack of polarized emission in the south- 
ern part suggests stronger Faraday depolarization 
in the southern than in the northern part of M33 
(see Sect. 3.6). In the northern arms, the degree of 
polarization P is larger than 20%, with a maximum 
of 30% on the arm IV N and also in a region close 
to NGC604 in IN. 



3.4. Distribution of total spectral index 

Because of their large frequency range, the total 
intensities at 20 cm and 3.6 cm lead to a reliable 
spectral index map. After smoothing to the angu- 
lar resolution of 90", the total intensity maps at 
20 cm and 3.6 cm were restricted to a common size 
and interpolated onto the same grid. The spectral 
index, a (I^j^~"), was only computed for pixels 
with a flux density of at least three times the rms 
noise a at both frequencies. In the resultant spec- 
tral index map (Fig. 9), the spiral arms and star 
forming regions are well visible as flat spectrum re- 
gions. Steeper spectra occur in the outer regions 
and in between arms. Fig. 10 shows that the total 
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Table 5. Integrated flux densities for B < 8.5 kpc (35'). 



A 


HPBW 


St 


Spi 


P 


Sp/(180") 


P(180") 


(cm) 


(arcsec^) 


(mJy) 


(mJy) 


% 


(mJy) 


% 


3.6 


120 


779 ± 66 


48 ±5 


6.1 ±0.8 


40 ±4 


5.1 ±0.7 


6.2 


180 


1284 ±135 


89±4 


6.9 ±0.8 


89±4 


6.9 ±0.8 


20 (combined) 


51 


2768 ± 63 


167 ±25 


6.0 ±0.9 


123 ±3 


4.4 ±0.1 



_ % \ * — \ r 



3 

a 41. 
5 




, #■ 



I r I I ^ 



Fig. 9. Total spectral index map of M33 as deter- 
mined between the two wavelengths 3.6 and 20 cm. 
The angular resolution is 90". For the total inten- 
sity maps, 3cr level was considered to obtain this 
spectral index map. 



spectral index has a mean value of 0.7. The error in 
the total spectral index varies from 0.01 in the star 
forming regions to 0.23 in the outer parts. 

The HII complexes IC 133, NGC604, and 
NGC 595 have spectral indices < 0.1 and even neg- 
ative values, indic ating optically thick c onditions in 
these HII regions. Johnson et al.l ( 2001 ) found sev- 
eral optically thick HII regions (at 6 cm and 20 cm) 
with the two brightest ones in IC 133. 

The most inverted spectrum is due to a source 
at RA = 01^ 32"^ 42.06±0.14" and DEC = ±30° 54' 
24.8±1.8" which has a spectral index of -0.6±0.1 
(l3.6cm = 3.6 ±0.2 and Isocm = 1.24 ±0.07). This 
source, which does not exist in the catalogues 
of radio sources and HII regions of M33 (e.g. 



m ICordon et al]|l999l ). could be a variable radio 
source or a giant optically thick HII region behind 
a high concentration of dust. However, the latter 
possibility is unlikely since no IR emission is ob- 
served from that region. 

We obtained the radial distribution of the spec- 
tral index a using the integrated total flux densities 
calculated in rings with a width of 0.5 kpc in the 
galactic plane. Fig. 11 presents the radial profile of 
the spectral index obtained from pairs of the inte- 
grated flux densities at 3.6 cm, 6.2 cm, and 20 cm 
considering one a level. There is a general increase 
in a from the ce nter to a radius of about 7 kpc in 
all three profiles. iBuczilowskil (|l988l ) explained this 
increase as a general decrease of the thermal frac- 
tion with radius. However, with our hi gh resolution 
study we detect structures not seen by iBuczilowski 
(|1988l ) : a prominant bump around i? = 3 kpc and 
a steeper increase in a beyond i? = 4kpc. This in- 
dicates that the thermal fraction decreases faster 
in regions at i?>4kpc than at smaller radii. The 
bump is located in the ring 2.5 < i? < 3 kpc domi- 
nated by interarm regions which are expected to be 
weak in thermal emission. 

Assuming that the radial increase of the spec- 
tral index is related to a decrease in the thermal 
fraction, and considering that the thermal emis- 
sion contributes more to the radio continuum emis- 
sion at smaller wavelengths, one expects more vari- 
ations (or faster radial increase) of the spectral in- 
dex obtained at small wavelengths. This is indeed 
visible in Fig. 11: a(3.6,6.2) increases faster than 
a(6.2, 20) and a(3.6, 20) at i? > 4 kpc. It is also seen 
that the larger the wavelengths, the smaller the 
spectral index (or flatter energy spectrum at larger 
wavelengths) . This is surprising as one expects that 
the spectrum of the radio continuum emission be- 
comes steeper at larger wavelengths because of the 
synchrotron emission. This could hint to the non- 
thermal energy losses at larger wavelengths. The 
reason may be the same as of the significant flatten- 
ing of the integrated spectrum of M33 at frequen- 
cies lower than 900 MHz : free-free absorption of the 
non- thermal emission by a cool (< 1000 K) ionized 
gas ( Israel et al. 19921 ). The presence of the opti- 
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mean = 0.7 
median = 0.70 : 
nms = 0.23 




Fig. 
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total spectral index 

10. Histogram of the spectral index distribu- 
determined from total intensities 



Table 6. Integrated total flux densities from 
Effelsberg sm-veys {R < 50'). 



A (cm) 


u (MHz) 


S (mJy) 


35.6 


842 


5377 ±1217' 


21.1 


1420 


2990 ± 440^ 


20 (21) 


1425 


2768 ± 63^ 


17.4 


1720 


2714 ±254' 


11.1 


2702 


1683 ±168' 


6.3 


4750 


1100 ±167' 


6.2 


4850 


1300 ± 135^ 


3.6 


8350 


779 ± 66^ 


2.8 


10700 


482 ± 150' 



' Buczilowski 
^ this paper 



T ^ - )Ka(3,6 
2-5: oa(3,20 
Aa(6,20 

2.0 - 
1 .5 7 
1 .0 - 



.5 



0.0 



2 4 6 8 

R (kpc) 

Fig. 11. Total spectral index averaged in rings of 
0.5 kpc width in the galactic plane versus the ra- 
dial distance from the center of M33. For the total 
intensity maps, one a level was considered before 
integration. 



cally thick HII regions ( Johnson et al. 200 ll ) could 
also enhance this trend. 

Away from star forming regions and in the outer 
parts of the galaxy, the thermal fraction is neg- 
ligible and the observed spectral index a can be 
taken as the non-thermal spectral index, a„. From 
Fig. 9, we obtain a„ ~ 1.0 ± 0.1, which is equa l 
to that found in M31 bv lBerkhuiisen et al.l (|2003l ). 
Assuming that the non-thermal spectral index is 
constant across MBS, we find average thermal frac- 
tions of 0.52 ± 0.15, 0.46 ± 0.12, and 0.18 ± 0.06 
at 3.6 cm, 6.2 cm, ai id 20 cm, respecti vely (follow- 
ing the method by iKlein et aL 1984 ). How these 
values change if the non-thermal spectral index is 
not assumed to be constant across the galaxy will 
be discussed in paper II (Tabatabaei et al. in prep.). 



V (MHz) 

Fig. 12. Radio spectrum of the integrated emission 
from M33 using the data listed in Table 6. 

3.5. Integrated radio continuum spectrum 

The integrated flux densities for i? < 50' (12 kpc) 
from Effelsberg surveys at 35.6, 21.1, 17.4, 1 1 .1, 6.3 , 
and 2.8cm were presented by iBuczilowskil ()l988l ). 
They obtained a spectrum with a spectral index of 
0.91 ± 0.13, steeper than foimd in many other spiral 
galaxies. 

We rederived the integrated spectrum of M33 
using three more data points from the new sur- 
veys. We calculated the integrated flux densities 
for the 3.6, 6.2, and 20 cm maps for the same 
area of i?<50'. The 3.6 cm and 2.8 cm maps were 
integrated up to i? = 44' and R = 22.5', respec- 
tively, providing a lower limit to the total flux den- 
sity. However, extrapolated values up to R = 50' 
show that the missing flux densities a re within er- 
rors. O ur values together with those of [Buczilowski 
" l988f) (see Table 6) lead to an integrated spectrum 



(Fig. 11) with a spectral index of 0.72 ± 0.04, which 
is in agreement with the mean spectral index ob- 
tained in Sect. 3.4 using a different method. 
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Table 8. Northern and southern integrated flux densities in total and polarized intensity plus degree of 
polarization. 



A 


qN 

iJ'J' 






qN 
Opj 


Opj 






pi 




cm 


(mJy) 


(mJy) 




(mJy) 


(mJy) 




% 


% 




3.6 


412 ± 18 


366 ±21 


1.12 ±0.08 


20 ±2 


20 ±2 


1.00 ±0.14 


4.8 ±0.4 


5.5 ±0.6 


0.9±0.1 


6.2 


631 ± 63 


653 ± 70 


0.97 ±0.14 


53±4 


36 ±2 


1.47±0.14 


8.4 ±1.0 


5.5 ±0.6 


1.5±0.2 


20 


1268 ± 73 


1435 ± 83 


0.88 ±0.07 


81 ±5 


42 ±3 


1.93 ±0.18 


6.4±1.1 


2.9 ±0.3 


2.2 ±0.4 



Table 7. Exponential scale lengths. 



A 



/ii<4kpc (kpc) /p>4kpc (kpc) 



20 cm 
6.2 cm 
3.6 cm 
24 /im 
70 fim 
160 



5.44 ±0.27 
3.73 ±0.21 
3.48 ±0.53 
2.07 ±0.23 
2.04 ±0.19 
2.51 ±0.19 



2.74 ±0.10 
2.35 ±0.14 
1.29 ±0.03 
1.35 ±0.09 
1.85 ±0.16 
1.63 ±0.09 




R (Kpc) 



Fig. 13. Radial profiles of total intensity (surface 
brightness) at radio and IR wavelengths. 

3.6. Exponential scale lengths 

Fig. 13 shows the radial distributions of the total 
radio and IR (Spitzer 24, 70, and 160 /im) intensities 
in M33 with galactocentric radius. There is a break 
in the profiles near i? = 4kpc, which is best visible 
at 3.6 cm and 24 /im. All the profiles can be well 
fitted by two exponential functions of the form 

I(R) = / ^oexp{-B./lR<4kpc) 
\ l4exp{-R/lR>4kpc) 

where, /q and are the intensity at i? = and 
R = 4 kpc, respectively, and I is the exponential 
scale length. 

We derived the exponential scale lengths of the 
radio emission at the three wavelengths along with 
those of the IR emission (Table 7). Generally, the 
scale lengths are larger in regions at i? < 4 kpc 



i? < 4 kpc 
i? > 4 kpc. 



(1) 



than at i?>4kpc. This indicates that the surface 
brightness decreases faster with increasing radius at 
R > 4 kpc at all wavelengths. Sinc e the 24 /xm emis 



sion (and also 3.6 cm emission, iTabatabaei et al 



20071) traces the star forming regions, this con- 
firms that the main star forming regions in M33 
are spread over an area with a radius of i? < 4 kpc 
without a dominant 'nuclear concentr ation' (differ- 



ent f rom the situation in NGC6946, iWalsh et al 
2002h . 



Table 7 also shows that the scale length in- 
creases with increasing wavelength in the radio 
regime. This can be explained by the fact that at 
longer wavelengths the non-thermal emission be- 
comes more dominant, and the influence of radial 
diffusion of cosmic ray electrons becomes more im- 
portant, leading to larger scale lengths. 

3.7. North-south asymmetry 

Fig. 7 shows that the polarized emission at 20 cm 
is much weaker in the southern than in the north- 
ern half. To investigate this asymmetry, we calcu- 
lated the integrated total flux density in both to- 
tal and polarized intensities separately over each 
half (Table 8). The average degrees of polarization 
are also given in Table 8. At all three radio wave- 
lengths, the beam width is 180". The ratio of the 
northern to the southern integrated total flux den- 
sity in polarized intensity, Spj/Spj, and also the 
ratio of the average degree of polarization, P^/ P'^, 
increases with wavelength. /S^, does not change 
from 3.6 cm to 6.2 cm and from 6.2 cm to 20 cm 
within the errors, but, it decreases slightly from 
3.6 cm to 20 cm. A wavelength-dependent asymme- 
try in the polarized intensity and average degree of 
polarization was also found in other galaxies e.g. 
NGC6946 (Beck, submitted to A&A). The inter- 
pretation is that there is an asymmetry in Faraday 
depolarization along the line of sight that may be 
caused in several ways: 

- asymmetry in the distribution of the thermal emis- 
sion, 

- asymmetry in the uniformity of the magnetic field 
strength, 

- special large-scale structure of the regular mag- 
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netic field, like e.g. a large-sca le axisymmetric heli- 
cal magnetic field fjUrbanik et al. ,1997 ). 
We will investigate these possibilities after the sep- 
aration of thermal/non-thermal components in a 
forthcoming paper (Tabatabaei et al. in prep.). 

4. Discussion and conclusions 

Radio observations with high resolution and sen- 
sitivity enabled us to obtain much more informa- 
tion about the distribution of the linearly polar- 
ized emission in M33 than in previous studies, es- 
pecially, at high frequencies where the Faraday ef- 
fects are negligible. Using the new data at 3.6 cm, 
we found that the magnetic field is well ordered 
in a spiral structure with the general orientation 
of the optical spiral arms. No clear ordered mag- 
netic field-optical arm anti-correla tion, as observed 
in ot her gal axies like NGC6946 ( Beck fc HoernesI 



19961 ). M81 (iKrause et al 



IC342 (|Krause et al.lll989bl : lKrausdll993h . is found 



1989aD. and in parts of 



in M33. Detection of high degrees of polarization 
(larger than 20%), not only in the northern arms 
but also in the central extended region and in parts 
of the southern arms at 3.6 cm and 6.2 cm (for the 
first time), is another significant result of this study. 
Generally, we found a patchy distribution of the de- 
gree of polarization. The degree of polarization in 
the direction of the HII complexes is less than 5%. 
We also found a weak degree of polarization and 
some distortion in the regular magnetic field in the 
central IS, which could be linked to the high star 
forming activity in this region. 

We derived a spectral index map of M33 with a 
high angular resolution (90") (t he previous spectra l 
index map was derived at 462", iBuczilowsk il ll988D . 
Both the mean total spectral index (0.72 ± 0.07, ob- 
tained from the 3.6 and 20 cm maps) and the total 
spectral index obtained from the integrated spec- 
trum (0.72 ± 0.04) are consistent with t he result of a 
statist ical study on 56 spiral galaxies by Gioia et "ahl 
(|l982l ). 

The radial profile of the spectral index (Fig. 11) 
shows not only a general increase but also fluc- 
tuations at 2<i?<4kpc, confirming that the de- 
crease in the thermal fraction causes an increase 
in the spectral index. Considering Fig. 13, at 
2 < i? < 4kpc, we found opposite fluctuations in to- 
tal intensity and in spectral index. These fluctua- 
tions (in Fig. 13) are again stronger at wavelengths 
at which either the thermal emission is stronger 
(i.e. 3.6 cm) or the emission correl ates well with the 



thermal emission (i.e. 24 /im, see iTabatabaei et al 



l2007l ). The area at 2 < i? < 3 kpc includes the region 
between the main arm I (IN, IS) and the central 



extended region of the galaxy, corresponding to the 
small maximum (bump) in the spectral index profile 
and a small minimum in the total intensity profile. 
The area at 3 < i? < 4 kpc includes the arm I N and 
NGC604, corresponding to the small minimum in 
Fig. 11 and the small maximum in Fig. 13. 

Like the star formation radial profile (24 /xm ra- 
dial profile, see Sect. 3.6), the mostly non-thermal 
20 cm radial profile is fitted by a larger exponen- 
tial scale length at i? < 4 kpc than at i? > 4 kpc. 
This is a surprising result, indicating that even 
the non-thermal emission (which is less confined 
to the star forming regions than the thermal emis- 
sion) mimics the distribution of its parents, the 
star forming regions. A spatial coupling between 
the non-thermal emission and star forming re- 
gions was also found in other late-ty pe galaxies 
(|Dumke et al.ll2000l : [ciwzv et al.ll2007[) . Therefore, 
it would be highly interesting to derive the distribu- 
tion of the non-thermal emission for a more realis- 
tic non-constant spectral index of the non-thermal 
emission (Tabatabaei et al. paper II). 
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